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HELICOPTER ROTOR ROTATIONAL NOISE PREDICTIONS BASED 
ON MEASURED HIGH-FREQUENCY BLADE LOADS 

Robert N. Hosier 

Langley Directorate, U.S. Army Air Mobility R&D Laboratory 

Ramani Ramakrishnan 
The George Washington University 
Joint Institute for Acoustics and Flight Sciences 

SUMMARY 

The results of tests, conducted at the Langley helicopter rotor test facility, which 
simultaneously measured up to 200 harmonics of the fluctuating aerodynamic blade surface 
pressures and the far-field radiated noise are reported. These tests were conducted on a 
full-scale nontranslating rotor system. The experimental apparatus and procedures, data 
reduction and analysis, generalized loading results, and rotational noise theory are dis- 
cussed. After their characteristics were determined, the measured blade surface pres- 
sures were converted to loading coefficients and used in an existing theory to predict the 
far-field rotational noise. The effects of the blade loading coefficients, chordwise loading 
distribution, blade loading phases, and observer azimuthal position on the noise calcula- 
tions are explained. A comparison of the calculated and measured noise shows generally 
good agreement through 300 to 600 Hz, depending on the discreteness of the loading spec- 
trum. These data form an important step in developing the ability to predict rotor noise 
from aircraft while in the design stage. 

INTRODUCTION 

The detectability and community acceptance of helicopters depend in a large part on 
the noise characteristics of the main and tail rotors. This rotor noise can be conveniently 
grouped into two general categories, broadband and discrete. Some of the sources of 
broadband noise are the tip and trailing-edge vortices and boundary-layer turbulence. 
Discrete noise consists of impulsive noise such as blade-vortex interaction and advancing- 
blade compressibility shock, thickness noise resulting from mass displacement near tran- 
sonic tip speeds, and rotational noise generated by thrust and drag forces. 

In most flight conditions the predominant rotor noise source is the rotational noise. 
For many years this noise was predicted from a knowledge of only the uniform aerodynamic 
loads (ref. 1) but, except at high tip speeds, agreement with measured acoustic data was 
poor beyond the first few harmonics. 


In the mid 1960 's theoreticians postulated that the rotational noise could be attributed 
to both the uniform and fluctuating aerodynamic lift and drag forces. Rotational- noise pre- 
diction theories, based on the fluctuating blade-force assumption, were then developed, 
and theoretical calculations of rotational noise were made (refs. 2 to 10). These theories, 
all of similar content, required a priori knowledge of the amplitude and distribution of the 
high-frequency fluctuating aerodynamic blade surface pressures to calculate the far-field 
rotational noise. However, no measurements of these pressures existed. (In refs. 6 
and 11 approximately 30 and 10 pressure harmonics, respectively, were measured, but for 
the purpose of the computations in this report on the order of 100 or more pressure har- 
monics are required.) Because of the lack of data, theoretical models of the fluctuating 
blade surface pressure spectra were developed and the agreement between calculated and 
measured rotational noise was found to be dependent on the loading model chosen. In ref- 
erence 6, five sound harmonics were computed and fell within 10 dB of measured data. 

In reference 12 the combined rotational and broadband noise from helicopter main 
and tail rotors is calculated from computed rotor performance parameters. The rota- 
tional noise calculations are made for hover and forward flight and the results agree with 
measured data within 10 dB up to about 10th sound harmonic. 

In reference 7 the rotational noise equations were developed in terms of blade loading 
coefficients (the ratio of a blade pressure harmonic to the uniform pressure) and 
chordwise-loading- distribution functions. Although spanwise loading distributions are 
also included, it is suggested that the loading can be concentrated at an effective blade 
radius. Reference 13 contains a critical review of the theories of references 6 and 7. 
Reference 14 suggests that the mechanisms for the fluctuating aerodynamic blade loads 
are inflow turbulence, gusts, rotor- fuselage interaction and rotor- wake interaction. 

This paper reports results from tests conducted at the Langley helicopter rotor test 
facility which simultaneously measured both the fluctuating-aerodynamic-blade surface 
pressures and the radiated rotor noise. In some conditions, as many as 200 pressure har- 
monics were measured and used in conjunction with the theory of references 7 and 8 to 
predict the far-field noise. Specific attention is given to the experimental equipment, 
instrumentation, and procedures, the data reduction and analysis, the generalized loading 
results, the rotational noise theory, a discussion of the blade pressure and loading data, 
and a comparison of measured and calculated rotational noise based on the measured fluc- 
tuating blade loading coefficients. The importance of chordwise loading distributions, 
blade-loading phase angles, and observer azimuthal position are also discussed. 

SYMBOLS 

The units for the physical quantities defined in this paper are given both in the Inter- 
national System of Units, SI, and the U.S. Customary Units. 
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A 

a 0 

B 

C 0 

c L 

Crp 

C 

D 

D F 

F 

Fi 

*0 

Jq(Z 0 ) 

J - ^ 

k d 

Kl 

L 

Lp 

Lt 


rotor disk.area, (ft^) 
speed of sound, m/sec (ft/sec) 
number of blades 

Fourier coefficient of zeroth order 
Fourier coefficient of m^ 1 order 
mean rotor lift coefficient, 6Crp/ cr 
rotor thrust coefficient, T^/pAV^ 
blade chord, m (ft) 
observer distance from source, m (ft) 
total rotor drag force, sin /3, N (lb) 
force vector, N (lb) 
i^h component of F 

force vector normal to y (see appendix A and fig. 20) 

Bessel function of the first kind of order q and argument Z 


nd F , 9 o 

drag operating constant , N/m z (lb/ft^) 

Ra Q M e 

lift operating constant ~ — sin a’, N/m^ (lb/ ft^) 

x\3-n 


load per blade per unit span, N/m (lb/ft) 


total rotor lift force L-p cos /3, N (lb) 


total rotor aerodynamic load, N (lb) 



effective Mach number, 0.8Mj- 


M e 

tip Mach number, fir^/aq 
m,mB sound pressure harmonic numbers 

N rotor shaft frequency, Hz 

P general blade pressure, N/m^ (lb/ft^) 

Pq,P s amplitude of qth or s*h fluctuating blade pressure harmonic, N/m^ (lb/ft^) 
P r ms root- mean- square amplitude of fluctuating blade pressure, N/m^ (lb/ft^) 

Pq amplitude of zeroth (uniform) blade pressure harmonic, N/m^ (lb/ft^) 

p acoustic pressure, N/m^ (lb/ft^) 

q mode number (q ± = mB ± s) 

R observer distance from rotor center, m (ft) 

r^ rotor tip radius, m (ft) 

r distance along rotor measured from rotor center, m (ft) 

r e effective rotor radius 0.8rt, m (ft) 

SP sound pressure, N/m^ (lb/ft^) 

SPL sound-pressure level, dB (re 0.00002 N/m^) 

SP m B amplitude of mBth sound pressure harmonic, N/m^ (lb/ft^) 

Sq jump function (Sq = 0 for q < mB; Sq = 1 for mB sq) 

s pressure or loading harmonic number 

Tq, total rotor thrust, N (lb) 
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t 


time variable, sec 


V t rotor tip speed ftr^., m/sec (ft/sec) 

x observer coordinate vector, m (ft) 


XpXgjXg rectangular coordinates 

y source coordinate vector, m (ft) 



y limit and directivity function, mBJ q (ZQ) 

4 

A increment 


6( ) Dirac delta function 

7),% arbitrary phase constants, rad 

9 blade- observer azimuthal angle, rad 

X mB wavelength of mB th acoustic harmonic 

p density of air, kg/m 3 (slugs/ft 3 ) 

cr rotor solidity, 

nr t 

a' observer elevation angle, rad 


t time period, sec 

(pq,<p s phase angles of the blade pressure or blade loading harmonic, rad 
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X m B chordwise loading distribution spectrum function 

ip rotor azimuthal angle, rad 

angular frequency 27 tN 

Subscripts: 

diff differential 

g gage 

int integrated 

rms root mean square 

EXPERIMENTAL EQUIPMENT, INSTRUMENTATION, AND PROCEDURES 

To provide the previously unavailable high frequency aerodynamic blade surface 
pressures necessary to validate the rotational noise theories, a full scale test was con- 
ducted to simultaneously measure the fluctuating blade surface pressures and the far -field 
radiated noise. This section outlines the experimental equipment, instrumentation, and 
procedures used in this test. 

Equipment 

Helicopter rotor test facility .- The fluctuating blade surface pressures were meas- 
ured at the Langley helicopter rotor test facility shown in figure 1. Although the rotor 
used had a diameter of 10.65 m (35 ft), this facility can be used to study rotors up to 
21.3 m (70 ft) in diameter. It has a maximum rpm, thrust, and torque capability of 
approximately 400 rpm, 66.723 kN (15 000 lb), and 24.4 kN-m (18 000 ft-lb), respec- 
tively. The facility is equipped with 60 slip rings to transmit the data from the rotating 
blades to the ground recording station. 

Rotor system .- The rotor system used for the tests was a four-bladed hingeless 
rotor. Table I provides a general description of the rotor. A more complete description 
can be found in references 15, 16, and 17. 

Microphone array .- The radiated far-field noise was measured with two microphones 
located 45.72 m (150 ft) from the rotor axis as shown in figure 2. The first microphone 
was 9.14 m (30 ft) above the rotor disk plane; the second microphone was in the rotor 
disk plane. 
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Figure 1.- Langley helicopter rotor test facility. 
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TABLE I.- ROTOR SYSTEM PARAMETERS 


Type rotor system Hingeless 

Diameter, m (ft) 10.67 (35) 

Number of blades Four 

Blade chord, m (in.) 0.33 (13) 

Blade weight per blade, N (lb) 382.5 (86) 

Airfoil section Modified NACA 0012 

Blade taper 0 

Total blade twist (root to tip), deg -5 

Disk area, m^ (ft^) 88.5 (962) 

Solidity 0.0818 

Forward blade sweep, deg 1.4 


Performance measurement system .- The rotor performance parameters - thrust, 
shaft torque, rpm, and geometric blade tip pitch - were measured by a system of strain- 
gage bridges and potentiometers. The output of these sensors was monitored with digital 
voltmeters and manually recorded in a test log. 


Instrumentation 

Two data acquisition systems, one for blade surface pressure measurements and one 
for acoustic measurements, were used for these tests. Each system, with its accompany- 
ing subsystems, is described in this section. 

Blade surface pressure measurements. - The high-frequency fluctuating blade pres- 
sures were measured using the system shown as a block diagram in figure 3. This system 
may be divided into three distinct sections: pressure transducers, signal conditioners, 

and data recording equipment. 

Pressure transducers: The high-frequency fluctuating- blade surface pressures 

were obtained from commercially available semiconductor strain-gage absolute-pressure 
transducers which were sealed at atmospheric pressure (see fig. 4). These paddle type 
transducers had a flat frequency response from dc to 5000 Hz with a resonance frequency 
of 25 000 Hz. When excited by 10 V dc, the transducer produced an output of about 1 mV 
per 1149 N/m^ gage (6 mV/psig). The entire transducer, including the bridge circuitry, 
was 1.27 cm (0.5 in.) long and 0.0889 cm (0.035 in.) thick with a diaphragm diameter of 
0.635 cm (0.25 in.). The transducers were mounted on the surface of one rotor blade as 
shown schematically in figure 5. Originally there were upper and lower surface trans- 
ducers at each station; however, the lower surface 30-percent- chord transducer failed 
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1 Pressure transducer 

2 Slip rings 

3 5-volt dc power supply 

4 Oscillator 

5 Balance box 

6 Filter system 

7 Amplifier 

8 Signal monitor 

9 Tape recorder 

10 Oscilloscope 

11 Galvanometer driver 

12 Oscillograph 

13 Time code generator 

14 Blade position indicator 



Figure 3.- Block diagram of blade pressure data acquisition system. 
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L-72-1546.1 


Figure 4.- Photograph of surface of pressure transducer used for 
blade pressure measurements. 
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(b) Schematic of transducer mounting. 

Figure 5.- Schematics of location and mounting of surface 
pressure transducers. 



before data acquisition began. To relieve case strains which may have developed because 
of blade bending, care was taken to achieve a cantilevered mounting by not bonding the 
transducer under the diaphragm area. 

Each pressure transducer was calibrated with a commercially available discrete 
frequency acoustic calibrator at a level of 62 N/m^ gage (0.009 psig) rms and at a fre- 
quency of 250 Hz. The calibrations were performed twice daily several days prior to each 
test. The mean level of these calibrations was taken to be the nominal calibration level. 

In order to achieve a constant amplitude calibration, an audio oscillator was used to elec- 
tronically simulate the mean level of the dynamic calibrations. This electrical calibration 
was introduced at the input to the balance box and recorded on FM magnetic tape as the 
transducer amplitude calibration. 

Signal conditioners: The output of the pressure transducers was passed through 

slip rings to a balance box where offsets due to changes in temperature, atmospheric pres- 
sure, and variation of individual gage zero balance were nulled. The slip ring noise, 
always of concern in this type of experiment, was found to have a maximum overall rms 
level of 50 pV or 57.45 N/m^ gage (0.0083 psig) based on the sensitivity of the pressure 
transducers. A typical time history and spectrum of the slip ring noise are shown in fig- 
ure 6. From this figure, it is seen that the slip ring noise is concentrated below 60 Hz; 
but this is also the range where the aerodynamic pressures are of largest amplitude. 

The balanced output was then amplified to obtain sufficient levels for FM recording 
on magnetic tape. Two different filters were also applied at this point. The first was a 
10 000 Hz low- pass filter and the second was a dc blocking band- pass filter with a 
10 000 Hz cutoff. The low-pass filter was used to measure the dc (uniform) pressures 
while the band- pass filter was used to measure the dynamic pressures. Both filters had 
essentially the same dynamic response, but the latter allowed amplification of the dynamic 
data without exceeding the FM recording limit. The frequency response of the blade sur- 
face pressure measuring system was found to be limited by the response of the trans- 
ducers. Figure 7 shows the frequency response of the dynamic pressure measurement 
system in the dc blocking configuration. The 6 dB down point is at 10 000 Hz. Because 
of the transducer frequency response limit, all of the pressure data analysis was confined 
to 5000 Hz. 

Recording equipment: All blade pressure data were multichannel, FM recorded on 

1-inch magnetic tape at a speed of 152.4 cm/sec (60 ips). The carrier frequency was 
108 000 Hz with a recording bandwidth of dc to 20 000 Hz. 

Time-code generator: A standard IRIG A 10 000 Hz time code was recorded on a 
separate recorder channel. This time code provided time markers every 0.1 sec which 
were used to identify tape segments for digitizing the data. 
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Position indicator: The azimuthal position of the instrumented blade was recorded 

four times per revolution by a blade pulse position indicator. The pulse was obtained 
from a voltage applied to a reed switch and capacitor circuit triggered every 90° by a 
magnet attached to the rotor shaft. 

Acoustic measurements .- Figure 8 shows a block diagram of acoustic data acquisi- 
tion system. The microphones were commercially available, piezoelectric, ceramic type, 
having a frequency response flat within 3 dB over the frequency range of 12 Hz to 
10 000 Hz. The signal outputs from both microphones were multichannel, FM recorded 
on magnetic tape. The frequency response of the complete measuring system was flat to 
within 3 dB from 12 Hz to 10 000 Hz. Figure 9 shows a typical frequency -response curve 
for the measuring system. The entire sound measurement system was calibrated in the 
field prior to and after each test. Calibrations were performed with a commercially 
available discrete frequency calibrator with a 1000 Hz sine- wave signal at a sound pres- 
sure level of 124 dB. 



1 Microphone 

2 Preamplifier 

3 Booster amplifier 

4 Sound level meter 

5 10 000 Hz time code 

6 Voice microphone 

7 Blade azimuthal position indicator 

8 Oscillograph 

9 Oscilloscope 

10 dc amplifier 

11 Real time analyzer 

12 Tape-recorder monitor switch 

13 Tape recorder 

14 110- volt power supply 


Figure 8.- Block diagram of acoustic data acquisition system. 
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Frequency, Hz 

Figure 9.- Typical frequency- response curve for the entire acoustic data system. 

Experimental Procedure 

A regular check on all the measuring systems was performed before operation of the 
facility. In order to insure low ambient noise levels, the experiments were conducted 
during the predawn hours. Wind velocities during these tests did not exceed 2.24 m/sec 
(5 mph). 

Table II lists the test conditions studied. Recording of the data was begun only after 
insuring that steady state was reached for each test condition. Each filter condition was 
used with two sets of amplification, the first with the maximum possible gain and the sec- 
ond with the next lowest gain setting. The latter case served as a precautionary measure 
in the event that the maximum gain case caused the data to exceed the voltage limit of the 
FM recording channel. The maximum allowable amplification for each test case was 
determined from oscilloscope monitors. In addition, a real-time narrow band analyzer 
was used to provide on-line displays of the pressure and acoustic spectra. After assuring 
that systems were operating properly, approximately 45 sec of blade pressure and acous- 
tic data were recorded simultaneously for each test case. 
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TABLE II.- TEST CONDITIONS 


Condition 

rpm 

Blade passage 
frequency, 
Hz 

Mach number, Mg, 
at 80- percent span 

Lift, 

l F 

Torque 

D F r e 

Angle of attack, /3 , 
at 80-percent span, 
deg 

N 

lb 

N-m 

lb- ft 

1 

225 

15.0 

0.30 

14 612 

3 285 

6 345 

4 680 

10.23 

2 

225 

15.0 

.30 

22 743 

5 113 

15 253 

11 250 

15.15 

3 

261 

17.4 

.35 

31 714 

7 130 

22 300 

16 450 

15.40 

4 

263 

17.5 

.36 

6 792 

1 527 

3 661 

2 700 

4.78 

5 

265 

17.7 

.36 

12 530 

2 817 

5 626 

4 150 

7.40 

6 

355 

23.7 

.48 

18 553 

4 171 

7 864 

5 800 

6.45 

7 

355 

23.7 

.48 

44 613 

10 030 

23 726 

17 500 

12.03 


DATA REDUCTION AND ANALYSIS 

Oscillograph records of all blade pressure and acoustic channels were reviewed fol- 
lowing recording of the test data; data which exceeded the recording limit were eliminated 
from further consideration. Test conditions which provided the most intense rotational 
noise were chosen for further analysis from the remaining data. A 20- second period of 
data from each of these conditions was chosen for analog- to- digital conversion. Initially 
the data were digitized at a rate in excess of 10 000 samples per second to provide blade 
loading coefficients to 5000 Hz. Examination of these data indicated that because of the 
calm inflow environment (hover with a 2.24 m/sec (5 mph) wind) no discrete blade pres- 
sure harmonics were visible beyond 1000 Hz. More severe operating conditions would 
have produced more intense high-frequency discrete blade pressures. To allow more 
efficient use of the computerized Fourier transform routines (ref. 18), the data were 
redigitized at a rate of 4000 samples per second. Antialiasing low-pass filters with a cut- 
off frequency of 2000 Hz were applied to the data before digitizing. 

The primary input parameters for the rotational- noise prediction computer pro- 
gram — fluctuating blade loading spectra - were obtained by first digitizing the blade 
pressure data and then using a Fourier transform computer program (ref. 18) to obtain 
pressure spectra. These complex spectra, produced for each of the seven absolute pres- 
sure transducers (fig. 5), were then summed for each pair of upper and lower surface 
transducers to give a differential pressure spectrum at each chordwise- station. After 
normalizing their amplitudes by the uniform pressure the differential spectra (now loading 
spectra) were then either used directly in the rotational- noise prediction computer pro- 
gram or they were chordwise integrated and supplied to the computer program for noise 
prediction. The integrated spectrum was obtained by employing a spline interpolation 
(ref. 19) routine and a combined Simpson-Newton 3/8 rule integration procedure (ref. 20). 
The uniform pressure was first obtained by measuring the dc shift of each transducer 
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from its zero position on the oscilloscope records. Later the methods of references 21 
to 24 were used because changes occurred in the transducer electronic balance. 

■ Although the prediction program required only the ‘“blade loading harmonics (periodic 
data), fluctuating blade pressures, by nature, contain random ^components. To increase 
the accuracy of the estimate of the blade load data, that is, to increase the statistical 
degrees of freedom, spectrum averaging was employed. Table III indicates the statistical 
degrees of freedom and the 90-percent confidence limits associated with each test case. 


TABLE III.- NINETY- PERCENT CONFIDENCE LIMITS OF DATA ANALYSIS BY FAST 

FOURIER TRANSFORM PROGRAM 


Condition 

rpm 

Analysis 

bandwidth 

Number of 
averages 

Degrees of 
freedom 

Error 
band, dB 

1 

225 

1.25 

25 

50 

1.6 to -1.8 

2 

225 

1.25 

25 

50 

1.6 to -1.3 

3 

261 

1.09 

21 

42 

1.4 to -1.8 

4 

263 

1.10 

22 

44 

1.4 to -1.8 

5 

.265 

1.10 

22 

44 

1.4 to -1.8 

6 

355 

1.48 

29 

58 

1.3 to -1.5 

7 

355 

1.48 

29 

58 

1.3 to -1.5 


The acoustic data, which were measured simultaneously with the fluctuating blade 
pressures, were also spectrum analyzed. Thus, it was possible to directly compare the 
predicted and measured rotational noise. The statistical degrees of freedom and the 
90-percent confidence limits for the measured acoustic spectra are also given in table HI. 

GENERALIZED LOADING RESULTS AND DISCUSSION 

Prior to calculating the rotational noise, the loading data from the seven test condi- 
tions listed in table il were studied to identify trends and to collapse the data where possi- 
ble. This section reports the results of that effort. 

Time Histories of Blade Pressure Data 

A total of seven transducers were used to measure the gage pressure on the upper 
and lower blade surfaces at four chordwise stations at 80-percent span. A pair of trans- 
ducers was used at the 15- , 50- , and 75-percent chord stations; but because the lower sur- 
face transducer failed, only a single transducer was used at the 30-percent chord station 
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(see fig. 5). Figure 10 shows typical pressure time histories from the upper and lower 
surface gages at 15- and 50-percent chord. These data (from condition 6, table n) show 
the fundamental 1/rev oscillation with higher frequency oscillations superimposed. The 
periodic nature of these data is indicative of steady but nonuniform blade pressures. It 
can be seen that the upper and lower surface pressure peak-to-peak amplitudes are nearly 
equal at each chord station. From figure 10 it can be seen that the percentage of high- 
frequency blade loads increases from the 15-percent to the 50-percent station; however, 
the peak-to-peak amplitude is less at 50-percent chord than at 15-percent chord. It can 




(d) 



Figure 10.- Typical gage pressure time histories from upper and lower surface 
transducers at (a) 15-percent chord, upper surface, (b) 15-percent chord, 
lower surface, (c) 50-percent chord, upper surface, and (d) 50-percent chord, 
lower surface. 


19 



also be seen that the upper and lower surface pressures at each station are approximately 
180° out of phase. 

Root- Mean- Square Levels of Blade Pressure Data 

Table IV (a) shows the rms gage pressures at each of the seven stations and the mean 
lift coefficients f° r all the test conditions. Again it is observed that except for sta- 

tions 4 and 5 the upper and lower surface rms pressures at each chord station are nearly 
equal with a general decrease in level from leading to trailing edge. These data are 
plotted in figure 11(a). Except for conditions 6 and 7, the rms levels remain nearly con- 
stant at each station for all C^. For comparison the maximum overall rms system 
noise level is also shown in figure 11(a). 


TABLE IV.- ROOT- MEAN- SQUARE GAGE PRESSURES AND 
MEAN LIFT COEFFICIENTS 


(a) Absolute values 


Condition 

c i 

Root- mean-square gage pressure, N/m 2 , at - 


Station 1 

Station 2 

Station 3 

Station 4 

Station 5 

Station 6 

Station 7 

1 

0.619 

903 

641 

268 

813 

220 

199 

186 

2 

.938 

827 

620 

227 

510 

220 

206 

275 

3 

.991 

923 

744 

275 

751 

241 

241 

337 

4 

.210 

717 

648 

289 

434 

172 

158 

151 

5 

.383 

779 

620 

268 

496 

172 

193 

165 

6 

.317 

1048 

889 

406 

586 

227 

227 

262 

7 

.759 

1399 

1165 

482 

854 

255 

372 

351 

Condition 

C T 

Root- mean- square gage pressure, psi, at - 

L 

Station 1 

Station 2 

Station 3 

Station 4 

.• * 

Station 5 

Station 6 

Station 7 

1 

0.619 

0.131 

0.093 

0.039 

0.118 

0.032 

0.029 

0.027 

2 

.938 

.120 

.090 

.033 

.074 

.032 

.030 

.040 

3 

.991 

.134 

.108 

.040 

.109 

.035 

.035 

.049 

4 

.210 

.104 

.094 

.042 

.063 

.025 

.023 

.022 

5 

.383 

.113 

.090 

.039 

.072 

.025 

.028 

.024 

6 

.317 

.152 

.129 

.059 

.085 

.033 

.033 

.033 

7 

.759 

.203 

.164 

.070 

.124 

.037 

.054 

.057 


(b) Percent of uniform load 


Condition 

C L 

Root 

-mean-square gage pressure, % of uniform load, at - 

Station 1 

Station 2 

Station 3 

Station 4 

Station 5 

Station 6 

Station 7 

1 

0.619 

12.7 

98.9 

5.4 

368.7 

4.5 

126.1 

15.9 

2 

.938 

8.5 

28.7 

3.5 

107.2 

5.7 

41.1 

17.7 

3 

.991 

' 9.0 

29.9 

4.1 

126.7 

6.0 

43.2 

20.8 

4 

.210 

16.8 

42.9 

8.8 

40.4 

8.3 

62.2 

20.7 

. 5 

.383 

14.0 

138.5 

6.6 

75.0 

6.8 

311.1 

17.6 

6 

.317 

20.4 

115.2 

10.7 

73.9 

9.5 

165.0 

30.2 

7 

.759 

16.5 

74.5 

8.4 

523.6 

7.3 

110.2 

25.6 
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rms gage pressure, N/m & (psig) 


(a) rms gage pressure for all stations. 

Figure 11.- Mean lift coefficient Cl plotted against rms gage pressures and 
rms loading coefficients o! rmS} g. (Test conditions refer to table II and 
station numbers refer to fig. 5.) 


In table IV(b) and figures 11(b) and 11(c) the rms gage blade loading coefficients, 

P rms g 

ffrmSjg = —p — are shown for all the test conditions. The upper surface coefficients 
1 0,g 


(fig. 11(b)) maintain nearly the same shape as the upper surface gage pressures (fig. 11(a)); 
but now, with the exception of station 1, the level of the coefficients increases going from 
the leading edge to the trailing edge. This is caused by the small normalizing uniform 
pressures, P n , at the trailing portion of the airfoil as compared with P n _ at the 
leading portion of the airfoil. This effect will be seen later in both the differential blade 
loading coefficients and in the rotational noise calculations. The lower surface blade load- 
ing coefficients (fig. 11(c)) behave much more erratically than the lower surface pressures 
(fig. 11(a)) or than the upper surface coefficients (fig. 11(b)). Their distorted shapes are 
caused by very small and erratic uniform blade pressures on the lower surface. 
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Blade Pressure Spectra 


Typical pressure amplitude spectra from data measured at stations 1, 2, 5, and 4 are 
shown in figures 12(a), 12(b), 12(c), .and 12(d), respectively. By comparing figures 12(a) and 
12(b) and figures 12(c) and 12(d), it can be seen that the upper and lower surface spectra are 
quite similar at each chord station. About 50 blade loading harmonics are evident in fig- 
ures 12(a) and 12(b) through about 300 Hz. Beyond this point a transition to a broadband 
spectrum occurs, although there seem to be intermittent areas of discrete data at these 
higher frequencies. The cause of these tones is not known. 

Although not obvious from figure 12, it will be shown in connection with figure 15 
that the blade pressures at the 50- and 75-percent chord stations generally tended to have 
less discrete content than did the pressures at the 15- and 30-percent chord stations. 

This increased randomness over the aft portion of the chord may be caused by increased 
boundary-layer turbulence in that area. In these tests, no blade pressure harmonics were 
ever observed above 1000 Hz. However, tests conducted with more severe inflow condi- 
tions may produce blade pressure harmonics up to or beyond 1000 Hz. 

System Noise 

In order to determine the pressure data system electrical noise, the pressure trans- 
ducers were replaced by dummy bridges while leaving the remainder of the signal condi- 
tioning intact. The rotor system was then operated over the range of conditions shown in 
table II. The pressure data system output included both the slip ring and signal condition- 
ing noise. Figure 13(a) to 13(d) show the typical pressure equivalent system electronic 
noise spectra for the data from stations 1, 2, 5, and 4. It can be seen in this figure that 
the system noise is concentrated below 60 Hz. The harmonic separation is equal to the 
shaft frequency and the harmonic levels are less than 11.6 N/m^ gage (0.00166 psig). 
Comparison of figures 12 and 13 shows that a minimum signal to noise ratio of about 10 dB 
(a factor of 3) is maintained for the 15-percent chord transducer at frequencies below 
60 Hz. Above 60 Hz the signal-to-noise ratio exceeds 20 dB. For the 50-percent chord 
transducer the signal-to-noise ratio is less than 10 dB at some of the harmonics below 
60 Hz. Above 60 Hz the signal-to-noise ratio exceeds 20 dB. The signal-to-noise ratio 
at the 30- and 75-percent chord stations was similar to that at the 15- and 50- percent 
chord stations, respectively. These signal-to-noise ratios were judged to be acceptable 
for these tests. 
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(b) 15-percent chord, lower surface. 


1.4 X 10 



(c) 50-percent chord, upper surface. 
Figure 12.- Continued. 
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(d) 50-percent chord, lower surface 
Figure 12.- Concluded. 
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(a) 15-percent chord, upper surface. 

Figure 13.- Typical pressure equivalent system electronic noise spectra 

for stations 1, 2, 5, and 4. 








Equivalent gage pressure, N/m^ (psig) 



Figure 13.- Concluded. 



Comparison of Measured Blade Pressures and Acoustic Spectra 

The measured far-field acoustic spectrum generated by the blade pressures shown 
in figures 12(a) to 12(d) is presented in figure 14. Comparison of the pressure and acous- 
tic spectral characteristics indicates some similarity between the spectra. The acoustic 
spectrum shows one acoustic harmonic for every four blade loading harmonics. This four 
to one correspondence is caused by the four blades on the rotor. The transition region 
where the discrete spectrum begins to become broadband in both the blade loading and 
acoustic spectra occurs between 300 Hz and 400 Hz. In these tests the discrete spectrum 
rarely extended beyond 500 Hz in either the pressure or acoustic data. 

The far-field acoustic spectrum amplitude decreases more slowly in amplitude than 
the. loading spectrum. This is caused by the fact that each sound harmonic is generated by 
a multiplicity of blade loading harmonics (see section entitled "Rotor Rotational Noise 
Theory"). 



Figure 14.- Measured acoustic spectrum at microphone position 1 
(see fig. 2). Thrust, 18 554 N (41711b); C L = 0.317; 
rpm, 355. 


33 



Differential Blade Loading Spectra 


Differential loading spectra at each station were obtained by performing a complex 
addition of the upper and lower surface gage pressure spectra and then normalizing by the 
differential uniform pressure. (Note: Hereinafter, all blade loading coefficients and 

pressures are based on gage pressures; therefore, the subscript g will be dropped.) 

At the 30-percent chord station where only an upper surface pressure transducer was 
used, the differential coefficients were approximated by doubling the amplitude of the upper 
surface blade pressure. This is believed reasonable because the upper and lower surface 
pressures are comparable in magnitude, but the upper surface uniform pressure is much 
greater than the lower surface uniform pressure. Hence, this doubling process approxi- 
mates quite closely the realistic values of the blade pressure coefficients as^diff* 


In order to characterize the differential blade loading coefficients a? S)( jiff spectra 
obtained at each condition listed in table II were superimposed on the same plot for each 
chordwise station. The results are shown in figure 15. The bands from the 15- and 
30-percent chord stations are very similar in level and shape (compare figs. 15(a) and 
15(b)). At these two stations, the levels o g ranged from 65 percent at low fre- 
quencies ^at the 15 -percent chord station excluding 0 Hz where a g = lj to 
0.003 percent at high frequencies (the lowest level at the 30-percent station). The 
spectrum fall-off rate ranges from 7 to 10 dB per octave (beginning with the first har- 
monic) throughout the entire spectrum. This fall-off rate is considerably less than the 
15 dB per octave fall-off rate suggested in references 3 and 6 but agrees well with the 
7 dB per octave fall-off rate, over the first 30 harmonics, reported in reference 25. 

The maximum level of the first harmonic is 65 and 25 percent of the differential uni- 
form pressure Pq for the 15- and 30-percent chord stations, respectively. 

Figures 15(c) and 15(d), taken from 50- and 75-percent chord, respectively, show 
a more rapid transition to broadband loads and larger levels in the high-frequency por- 
tion of the spectrum. However, the level of the first harmonic remains nearly the same 
as at the forward chord stations. 


The reason for the increase in the high frequency spectrum levels at the trailing 
stations (when compared with the leading stations) can be understood from figure 16. This 
figure shows typical differential blade pressure (P Sj diff = <*s,diff x Pq diff) s P ec * :i ' a f rom 
15- and 75-percent chord. Although the levels of the low-frequency harmonics are less at 
75-percent chord than at 15-percent chord, the high-frequency broadband pressures are 
higher in level at the 7.5-percent chord position. Also, the Pq at 75-percent chord 
is 17 percent of the Pq at 15-percent chord. These combined effects result in much 
higher blade loading coefficients aig ^iff at the aft chord positions, especially at high fre- 
quencies. As will be shown later, the low fall-off rates, the quick transition to broadband 
loads, and the high levels at 50- and 75-percent chord have a direct effect in the rotational 
noise calculations. 
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(a) 15-percent chord. 

Figure 16.- Typical differential blade pressure spectra. 





Figure 16.- Concluded. 



Integrated Differential Blade Loading Coefficients 

Because of the variation of the differential blade loading coefficients from chord 
station to chord station, the following procedure was used to produce loading coefficient 
spectra representative of the entire 80-percent span station. The amplitude and phase of 
each differential blade pressure harmonic measured at each chord station were used in a 
spline interpolation routine (see ref. 19) to provide better definition of the chordwise pres- 
sure distribution. In performing the interpolation, the conventional assumption that the 
blade pressure amplitude is zero at the leading and trailing edges was made. The inter- 
polated blade pressure amplitudes and phases were then chordwise integrated using a 
Simpson-Newton 3/8 rule (see ref. 20). Integrated differential blade loading coefficients 
0 's were then formed by normalizing the amplitudes of the integrated pressure har- 

monics by the integrated differential .uniform lift Pq ^ff Figure 17(a) shows a typical 
full- chord integrated diff int spectrum. It is characterized by a predominance of low- 
frequency discrete tones which gradually make a transition into a continuous spectrum near 
200 Hz. The first harmonic level is 41 percent of the uniform load. The spectrum fall- off 
rate begins at about 12 dB per octave and then gradually flattens to a 6 dB octave rate. 

The flattening of the spectrum is attributed to the influence of the higher amplitude broad- 
band data at the 50- and 75-percent chord stations. 

In order to reduce the influence of the more broadband loading at 50- and 75- percent 
chord, the amplitudes of the blade pressures beyond 40-percent chord were set equal to 
zero after interpolation. The pressures were then integrated and normalized over the 
leading 40 percent of the rotor blade. Figure 17(b) shows a typical 40-percent chordwise- 
integrated blade loading coefficient spectrum. This spectrum is characterized by a pre- 
dominance of low-frequency discrete tones which gradually make a transition into a con- 
tinuous spectrum. However, the transition occurs near 300 Hz rather than 200 Hz as in 
the full- chord case. The first harmonic level is about 47 percent of the uniform load. 

This is comparable with the 41 percent observed in the full- chord case. The fall- off rate 
is about one- half of the full- chord case, beginning with a 6 dB per octave and increasing 
to 8 dB per octave. Occasional discrete spikes can also be observed in the continuous por- 
tion of the spectrum. These may be caused by periodic inflow turbulence, although there 
is no conclusive evidence to support this. Also, although the spectrum does tend to become 
continuous beyond 300 Hz, there are intermittent areas of discrete spectrum up to 700 Hz. 

Superimposed Integrated Differential Blade Loading Spectra 

The superimposed integrated differential a S; diff ,int spectra at 80 percent span 
from aU the conditions listed in table II are shown in figures 18(a) and 18(b). Figure 18(a) 
shows the superimposed full chord integration spectrum. It reflects the fall- off trend seen 
in the typical spectrum, a rapid initial fall-off with a rather rapid flattening of the fall-off. 
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(a) Full chord integration. 

Figure 17.- Typical integrated differential blade loading coefficients. 



100 



40 


(b) 40-percent chord integration. 
Figure 17.- Concluded. 





Frequency, Hz 

(b) 40-percent chord integration. 

Figure 18.- Superimposed integrated differential blade loading coefficients. 
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The superimposed 40-percent chord integrated differential a S) diff,int spectra are shown 
in figure 18(b). Again the fall- off trend observed in the typical spectrum is reflected in 
the superimposed spectrum. The spectrum begins with and sustains a 7 to 10 dB per 
octave fall- off. The fundamental blade-loading harmonic amplitude is seen to have a maxi- 
mum level of about 65 percent of the uniform loading in both the full and 40-percent inte- 
gration cases. 

ROTOR ROTATIONAL NOISE THEORY 
Background 

The classical theory on discrete noise radiation is that presented in reference 1. 
There the far-field acoustic analysis was based on the sound radiated from a uniform blade 
loading force under steady inflow conditions. The predictions resulting from reference 1 
agreed with measured rotational noise data only for rotors and propellers with high tip 
Mach numbers. At moderate and low tip Mach numbers the results diverged drastically 
from the measured data (see fig. 19), with the exception of the first harmonic. 

Reference 26 extended the work of reference 1 by incorporating forward- speed effects 
for a propeller. The theories of references 2 and 7 were the first to formulate the fluctu- 
ating aerodynamic blade loads concept for rotor rotational- noise prediction. In reference 6 
measured blade loading data were used to predict rotor rotational noise up to the 10th 
acoustic harmonic but good agreement was achieved for about the first 5 acoustic 
harmonics. 

In order to allow calculations of the rotational noise based on the blade loading data 

measured in this study, the theory of references 7 and 8 was used. This section outlines 

'\ 

the salient points of that theory. 

Fluctuating- Blade Loads Theory 

The objective of the analysis is to solve a wave equation of the form 

1 d 1 2 P 9 2 P aQ 9F i 92T ij 

a 2 3t^ ax? 91 9x^ 3x^3xj 

1 

The mass source monopole Q accounts for discrete thickness noise radiation; however, 
it is considered of secondary importance for low and medium tip speed rotors and will not 
be included in the analysis. The stress tensor Tij is a quadrupole source and, therefore, 
will also be neglected in the analysis. This leaves only the force F (the bar indicates a 
vector quantity) on the right-hand side of equation (1). The forces and rotor coordinate 
system used are shown in figure 20. 
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Figure 19.- Comparison of measured rotational noise with calculated values 
obtained by using the procedure of reference 1. 
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In deriving equation (1), it was assumed (refs. 2 and 27) that the fluid was an inviscid 
and nonconducting medium. The force field component R should include all the loads on 
the rotor system. In the theory of reference 7 the forces are considered to be concentrated 
at an effective rotor radius of 0.8 span. Figure 21 conceptually illustrates the components 
of these aerodynamic blade pressures (or forces) at an arbitrary effective radius r . It is 
seen that the pressures are composed of uniform and fluctuating components. Further- 
more, the fluctuating component contains both periodic and nonperiodic pressures. Only 
the uniform and periodic components are required to calculate rotational noise and, as 
will be seen later, use of nonperiodic data causes the calculated noise levels to deviate 
from the measured noise. A detailed solution of equation (1) for the far-field radiated 
sound pressure is presented (following the analysis of ref. 7) in appendix A. Highlights 
of that solution are now discussed. 
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Figure 21.- Components of aerodynamic blade pressures. 


If loading phases and observer azimuth angle 9 are ignored, the contribution of the 
sth blade loading coefficient to the mB^ 1 far-field sound pressure harmonic is (see appen- 
dix A, equation A(15)) 


— mB ,s 2 X mB ( K L 





( 2 ) 


Blade loading coefficients .- For purposes of interpretation, equation (2) can be 
divided into three quantities. The first quantity is the contribution of the fluctuating blade 
loads given by the blade loading coefficient Og = P s /Pq. As discussed earlier, the coef- 
ficients Og may come from individual blade measurement stations or may be integrated 
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across the chord to provide a spectrum representative of a particular spanwise station. 
Both station and integrated blade-loading- coefficient spectra are used later in this report 
to calculate the rotational noise. 

Chordwise loading distribution .- The second quantity of equation (2) is the loading 
distribution function x m B* This function could describe both the spanwise and chordwise 
distribution of the uniform and fluctuating blade loads. However, since the spanwise dis- 
tribution was ignored by concentrating all the loading at 0.8 span, only the chordwise dis- 
tribution of the blade loads was included in this analysis. Figure 22 illustrates a typical 
chordwise loading distribution. In figure 22(a) the resultant uniform and fluctuating loads 
are shown at a given instant of time. The mean area under the distribution represents the 
total load per blade. The magnitude of this load is absorbed in the operating parameters 
Kl and Kq. The shape of the distribution is contained in the function x m B> determined 
from a Fourier decomposition of the distribution as a function of time. Figure 22(b) shows 
a time history of the loading distribution as seen at an arbitrary, but fixed, point on the 
rotor disk. It will be noted that the reaction of the air is a reflection of the blade pressure 
distribution about the leading edge of the airfoil. Also, the time history has a periodic 

nature with a period r = -i-. The duration of the distribution is equal to the blade chord 

BN 

divided by the rotor speed at the disk point. Therefore, x m B is a function of the acoustic 
harmonic number mB. 



(a) Distribution at instant of time.' 



(b) Time history, of loading distribution. 
Figure 22.- Typical chordwise loading distribution. 
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In actuality each blade loading harmonic may have a unique chordwise distribution 
making the distribution function dependent on both the acoustic and loading harmonic num- 
bers. For the purposes of this study both the uniform and fluctuating blade loads were 
assumed to be chordwise distributed in the same manner. Table V, after reference 8, 
shows the chordwise distributions and associated y-functions used in this study. Although 
not reported here, the exact a s distributions were measured when the data for this report 
were acquired. 


Operating constants. - The remaining terms in equation (2) form the third quantity in 

the equation. The term K T = L-r, sin o' is called the rotor lift constant and is a 

D RaQ r 

function of the shaft frequency N, the observer distance R and observer elevation angle 

TJ Dp 

a', the speed of sound a^, and the total rotor lift force Lp. The term K D = is 

called the rotor drag constant and is a function of the total rotor drag Dp, the observer 


TABLE V.- CHORDWISE LOADING DISTRIBUTIONS AND ASSOCIATED X" FUNCTIONS 
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distance R, the speed of sound ag, and the effective Mach number Mg (usually the 
Mach number at 80-percent span). The terms yq = mBJq (mBMg cos a') and 

yq + = mBJq + (mBMg cos a') are functions of the acoustic harmonic number mB and 
the Bessel function Jq of argument mBMg cos o'. The indices q^. = mB ± s. The 
combined effect of these terms is to provide the net lift and drag contributions and overall 
directivity according to the elevation angle a'. 

General Discussion of Theory 

In developing equation (2) the following assumptions were made in reference 7: 

(1) Equation (1) is valid only for a nontranslating system. 

(2) Far-field approximations are made. 

(3) All the rotor loads are concentrated at an effective radius of 0.8 span. 

(4) The rotor inflow is steady but nonuniform; therefore, each blade experiences the 
same load fluctuations every revolution. (Although the inflow is steady and nonuniform, 
the rotor loads can be Fourier decomposed into steady uniform, or dc, loads and steady 
nonuniform, or fluctuating, loads.) See figure 21. 

(5) The loading distribution across the chord is described by generalized loading pro- 
files. (The Fourier transform of these loading distributions yields the x _ functions dis- 
cussed earlier in this section.) Typical distributions are: point loading, rectangular load- 
ing, half- cosine loading, triangular loading, and sawtooth loading. 

The first three assumptions are implicit in the theory but the last two are subject to 
experimental verification. The pressure- time histories presented earlier (see fig. 10) 
seem to indicate a steady but nonuniform inflow. Even so, for the purposes of this study, 
the blade loading spectra were made steady by the averaging techniques discussed in the 
section "Data Reduction and Analysis." This is seen in the periodic nature of the pressure- 
time histories and in the discrete nature of the spectra (see fig. 12). In the next section 
the effect of the various chordwise loading distributions on the radiated noise is discussed, 
and experimentally determined loading coefficients are used to compute the radiated rota- 
tional noise. 

From equation (2) it is seen that for a particular mB and s number the resulting 
sound pressure is generated by two modes, q_ and q + . These modes are defined as 
rotating sinusoidal pressure patterns with spinning speeds of mB and mB for 
the q_ and q + modes, respectively. The spinning effect causes each mode to move 
periodically toward and away from the observer, allowing each blade loading harmonic to 
contribute to the level of more than one sound harmonic. Viewed from the observer' s 
frame, each sound harmonic is composed of the contributions of many (in theory an infinite 
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number) loading harmonics with the strongest contribution coming from the loading harmonic 
whose order s is equal to the order of the sound harmonic mB. 


The final result from this multimodal analysis, including all loading phases 


the observer azimuthal angle d, and the retarded time 
eq. A(21)) 



is (see appendix A, 


0q±mB> 



+ -T SK L 5 'O ex « , (-i*mB) ( -i> 



For a particular sound harmonic mB, the Bessel function J q varies rapidly with 
loading harmonic number s. Figure 23 shows the result obtained in reference 27 where 
the magnitude of the Bessel function Jq(Zg) * 10"^ or Jq(Zg) =10“^ for specific 
values of the order and argument. For this work the criterion 



(where Zq = mBM e cos o’ and £ 3 indicates the integer part of 10 + 5/4 Zq) was 
chosen as the maximum Bessel function order because any Jq(Zg) with q > |10 + 5/4 Z^j 
were assumed to equal zero. Thus, the infinite summation in equation (3) can be replaced 
by a finite upper limit given by equation (4). Since q = mB ± s, the finite limit is inter- 
preted as defining a band of loading harmonics s, centered about s = mB (that is, 
q = 0), which contribute to each sound harmonic mB (see refs. 7 and 8 and appendix A). 

In reference 8 the modes q + = mB + s were dropped from the analysis. However, 
in reference 28 it was commented that the terms q + have a nonnegligible effect on the 
calculations at low acoustic harmonic mB. For this reason the modes q ' were retained 
in the present analysis. 
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Figure 23.- The effect of argument Zq and order q 
on the Bessel function Jq^Zp^. 


If, in equation (3), only the uniform loading (s = 0) is considered and all the phases 
are neglected (as indicated by the bar), equation (3) reduces to 

< 5 > 

This is the same result obtained in reference 1. 

COMPARISON OF CALCULATED AND MEASURED ROTATIONAL NOISE 

This section presents a comparison of rotational noise calculations, using the theory 
of references 7 and 8 and the measured blade loading data discussed earlier, with rotational 
noise measured at microphone positions 1 and 2 (see fig. 2). Data are shown for all seven 
conditions listed in table II. Before discussing the final results, however, data are pre- 


50 



sented which show the effect of various a spectra, chordwise distribution (x) functions, 
and blade-loading phase angle/observer angle on rotational noise calculations. 

Effect of Blade Loading Spectra on Calculated Rotational Noise 

In order to implement the theory of references 7 and 8, a knowledge of the blade 
loading coefficients Og is required. These coefficients could be obtained from measure- 
ments at a single chordwise station or they could come from an integration of the station 
loading coefficients. In order to determine which loading coefficients allowed the most 
accurate calculations of the rotational noise, differential blade-loading coefficient spectra 
from 15- , 30- , 50- , and 75-percent chord stations and from full and 40- percent chordwise 
integrations were used with a point loading distribution to predict the rotational noise. 
Typical results from these calculations (based on condition 6, table II) are presented in 
figure 24. Here calculations using each of the above mentioned blade loading spectra are 
compared with measured rotational noise spectra. (The circles represent the noise level 
at the frequencies mBN, harmonics of the blade passage frequency. It should be noted 
that the noise at these frequencies may not be caused by rotational noise beyond the 
rotational-broadband transition range.) It can be seen from figure 24(a) that the calcula- 
tions based on loading data from the 30- , 50- , and 75-percent chord station and the full 
chord integration did not provide good agreement with the measured noise. Calculations 
based on loading data from 15-percent chord and a leading 40-percent chordwise integra- 
tion provided good agreement with the measured noise throughout all but the highest five 
acoustic harmonics. In a majority of the calculations the 40-percent chordwise loading 
integration provided the. best agreement with measured rotational noise. The poor agree- 
ment using the 50- and 75-percent chord station and full chord integration data can be 
understood in terms of the data shown in figures 15 to 17. In these figures the strong 
influence of broadband loading over the trailing 50 percent of the airfoil was noted. Since 
the theory of reference 7 is meant to calculate discrete rotor noise from discrete loading 
spectra, the broadband content of the 50- and 75-percent chord station and full chord inte- 
gration spectra causes the calculated rotational noise to not agree with measurements. 

Conversely, the good agreement between measurements and calculations using meas- 
ured blade loading coefficients from the 15-percent chord station, or integrated over the 
leading 40 percent of the airfoil, may be attributed to a more accurate representation of 
the discrete blade loading spectrum. There were occasions where calculations using the 
40-percent integration did not provide the best agreement with measurements. Fig- 
ure 24(b) shows such a case. Here, the loading data from 30-percent chord provided the 
best agreement with measured noise. However, the 40-percent chordwise integration pro- 
vided the next best agreement. Thus, because the leading 40-percent chordwise integrated, 
differential loading spectrum provided a good representation of the discrete blade loads and 
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(a) Microphone position 2. 

Figure 24.- The effect of various blade loading spectra on calculated rotational 
noise, obtained by using a point loading distribution. 


52 



90 




Measured 

40-percent chord integration 
Full chord integration 
15-percent chord station 
30-percent chord station 
50-percent chord station 
75-percent chord station 


20l— 


0 100 200 300 400 500 600 

Frequency, Hz 

I I J I L I L_ I I 

4 16 28 40 52 64 76 88 100 

mB 

(b) Microphone position 1. 

Figure 24.- Concluded. 
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because it seemed to provide good agreement between calculated and measured rotational 
noise, it was used to compute the rotational noise for all the conditions listed in table II. 
These data will be shown later. 


Effect of Chordwise Distribution of the Blade Loading Coefficients 
on Calculated Rotational Noise 

In discussing the theory of references 7 and 8, it was mentioned that the chordwise 
loading distribution function influences the rotational noise calculations. In order 

to determine the strength of this effect, the leading chordwise integrated, differential load- 
ing spectrum was distributed according to each of the functions shown in table V. Typical 
results (based on condition 6, table n) are shown in figures 25(a) and 25(b). At the low 

acoustic harmonics, all the distributions seem to provide the same calculated noise level. 

As the harmonic number increases, however, the influence of each distribution becomes 
apparent. In figure 25(a), typical of a majority of the cases, it was observed that the point 
loading distribution provided the best agreement between calculated and measured rotational 
noise. The rectangular distribution always gave the poorest agreement, while the remain- 
ing distributions occasionally provided the best agreement with measured data. One such 
case is shown in figure 25(b). Here, three distributions provide better agreement with the 
measured noise than did the point loading distribution. These data indicate that the loading 
distribution can affect the calculation results particularly at the higher acoustic harmonics. 


Effect of Blade- Loading Phase Angle and Observer Angle 
on Calculated Rotational Noise 

In the absence of any conclusive experimental evidence, previous experimenters have 
assumed either zero or random blade loading phases (refs. 7 and 8). In reference 29 it is 
concluded that neither the blade loading phases nor the observer angle 9 is important for 
calculating rotational noise for a hovering rotor. In the strictest sense the latter statement 
is certainly true if a hovering rotor is construed to have only steady, uniform blade loading. 
However, nonuniformity in the loading, even for the nontranslating rotor studied in this 
report, will cause a variation in the radiated rotational noise with observer angle 9. 
Theoretically, the effect of observer angle and blade loading phase (p can be seen in 
exponential terms of the general form expj|j^q0 + ^|q ±rn B|j] ( see e< L (3))* Here q is 
a blade-loading model index. This expression implies that the value of the loading phase, 
0|q ± mB| will affect the calculated acoustic harmonic SP m g in the same manner as the 
product q0. Figure 26 illustrates this effect conceptually. Here, two loading modes, one 
high order and one low order, are shown in planform superimposed on the rotor disk. The 
observer position and its associated angle are shown in reference to the line where the 
Fourier analysis of the blade loading began (appendix B explains the experimental proce- 
dure used to determine 6). All loading phases are referenced to this line. Considering 
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(a) Microphone position 2. 


Figure 25,- The effect of various chordwise loading distributions on calculated 
rotational noise, obtained by using a 40 percent chord integrated loading 
spectrum. 
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(b) Microphone position 1. 
Figure 25.- Concluded. 
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Figure 26.- Conceptual view of typical low-order and high-order blade loading 
modes and their orientation with respect to observer. 

the low- order mode, it can be seen that the noise radiated would be greatly affected by a 
phase change of, say, 7t/2. The acoustic radiation from the high-order mode, however, 
can be expected to be much more nondirectional than for the low- order mode. Therefore, 
a phase change of 7t/2 would be expected to have little effect on its acoustic radiation 
pattern. It seems, then, that the blade loading phase and observer angles are most impor- 
tant in determining the level of the low- order sound harmonics. 

In order to test this effect, the Fourier transform program described in reference 18 
was used to compute the phases of the measured blade loading harmonics. The typical 
result, as shown in figure 27, was that the phase spectrum was concentrated in a band of 
±25° about 0°. Four trial cases were then conducted for a typical test condition. Three 
of the four cases consisted of using arbitrary, but constant, phases for the 0's spectra; 
the fourth case used the actual measured a s phase spectrum. The results from these 
four cases are presented in figure 28. Except for the first acoustic harmonic, all the cal- 
culated harmonics are coincident. This effect could be expected to spread to a few higher 
harmonics if the rotational speed were lower. 

Final Calculation Results 

Because the previously cited rotational noise calculations, based on 40- percent chord- 
wise integrated, differential blade loading coefficients (a s in ^ and a point chordwise 
loading distribution ^X m B = 1) > seemed generally to provide the best agreement with 
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Loads. 





<) - Predicted with 85° constant phase 
□ - Predicted with 285° constant phase 
©- Predicted with 0°constant phase 



using a rectangular loading distribution. 



measured noise, similar calculations were performed for all the conditions listed in 
table II. The resulting calculated and measured spectra at microphone positions 1 and 
2 (see fig. 2) are displayed in figures 29(a) to 29(n). (In these figures the circles indicate 
the measured noise levels at frequencies mBN, harmonics of the blade passage frequen- 
cies. Therefore, at frequencies above the rotational- broadband transition range, they are 
not the rotational noise level.) Some general comments can be made regarding these 
results: 

(1) The use of the point loading distribution and integrated blade loading coefficients 
is not to be construed as a criterion for rotational noise prediction. Their use in these 
calculations was strictly an engineering judgment since they seem to provide the best 
agreement with measured rotational noise most of the time. 


(2) Generally, the calculated and measured rotational noise spectra show good agree- 
ment through 500 Hz. Beyond 500 Hz the calculated and measured spectra begin to diverge 
in most cases. This may be attributed to the fact that at these frequencies the compact 
source model begins to break down, that is, the criterion * mB/ /2 < c is approached 

= c at a frequency of 523 Hzj. It is at these frequencies that the shape of the chord- 
wise distribution becomes most important. In order to accurately calculate the rotational 


noise levels for 


l mB 


< c, it may be necessary to use a unique chordwise distribution for 


each blade loading coefficient ^iff Implementation of such a technique was beyond 
the scope of this work, although the information is implicit in the data. 


(3) There are several examples (figs. 29(a), 29(b), and 29(h)) where the calculated and 
measured spectra diverge much below 500 Hz. This is not caused by a deficiency in the 
theory but seems to be a direct result of an early transition of the loading spectrum from 
discrete to broadband content. 


(4) In a number of cases (figs. 29(a), 29(b), 29(g), 29(h), 29 (j), 29(1), and 29(m)) the 
first-order harmonic is either under or over predicted by 7 to 14 dB. This is believed to 
be caused by the combined effect of small errors in the measurement of the blade observer 
angle 9 and loading phases cp s , plus a breakdown in the far -field approximations. The 
effect of 8 and (p s on the first-order harmonic amplitudes was discussed earlier. 
Because the microphone positions were 45.72 m (150 ft) from the rotor axis, even at the 
highest RPM tested they were only slightly more than three acoustic wavelengths away.. 

This would put them on the boundary between the near field and far field. Thus, better 
agreement at the lowest acoustic harmonics may, in this case, require the use of near -field 
terms in the analysis. 

(5) The generally good agreement between the measured and calculated rotational 
noise, based on loading measured at 80-percent span, suggests that the aerodynamic activity 
in that area is representative of the entire rotor system, at least as far as rotational -noise 
generation is concerned. 
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(a) rpm = 225; L_ = 14 612 N (3285 lb); C T = 0.619; a' = 0° (microphone position 2). 
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(b) rpm - 225; Lj, - 14 612 N (3285 lb); C-^ = 0.619; o' = 11°' (microphone position 1). 

Figure 29.- Comparison of calculated rotational noise, obtained by using 40-percent chord- 
wise integrated blade loading coefficients and point loading distribution, with measured 
noise for all conditions listed in table n. 
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(e) rpm = 261; Lp = 31 714 N (7130 lb); C L = 0.991; a' = 0° (microphone position 2). 
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(f) rpm - 261; Lp = 31 714 N (7130 lb); Cp = 0.991; o' = 11° (microphone position 1). 

Figure 29.- Continued. 
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(g) rpm = 263; Lp = 6792 N (1527 lb); = 0.210; o' = 0° (microphone position 2). 
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(h) rpm = 263; L p = 6792 N (1527 lb); C L = 0.210; o' = 11° (microphone position 1). 

Figure 29.- Continued. 
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(i) rpm - 265; Lp - 12 530 N (2817 lb); = 0.383; a’ = 0° (microphone position 2). 
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(j) rpm = 265; Lp = 12 530 N (2817 lb); C L = 0.383; o' = 11° (microphone position 1) 

Figure 29,- Continued. 
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(k) rpm = 355; Lp = 18 553 N (4171 lb); Cp = 0.317; cr' = 0° (microphone position 2). 



(1) rpm - 355; Lp = 18 553 N (4171 lb); Cp = 0.317; cr' = 11° (microphone position 1). 

Figure 29.- Continued. 
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CONCLUDING REMARKS 


In tests conducted at the Langley helicopter rotor test facility, simultaneous meas- 
urements were made of the fluctuating aerodynamic blade surface pressures and the radi- 
ated rotor noise. The measured fluctuating surface pressures were then used in Wright's 
theory to predict the rotational noise. Comparison of the predicted and measured rota- 
tional noise showed that the accuracy of the predictions was directly dependent on two 
things: the number of discrete loading harmonics measured and the chordwise distribu- 

tion of these loads. In general, test conditions which generated a large number of loading 
harmonics and which were assumed distributed at a point provided good rotational- noise 
predictions through 500 to 600 Hz. Within 1 to 5 dB, this is evidence that the compact 
source assumptions were valid for the operating conditions studied. At higher tip speeds, 
more realistic loading distributions would probably have to be considered. 

Examination of the measured blade pressures revealed that the upper and lower sur- 
face fluctuating pressures were quite similar in spectral amplitude and shape. Compari- 
son of the spectra at different chordwise stations, however, showed that the forward chord 
position contained more pressure harmonics through higher frequencies (up to 500 Hz) 
than did spectra from aft chord positions. The overall rms pressures were also greater, 
nearly 1500 N/m^ gage (0.217 psig), at forward chord positions. The aft chord positions 
generally contained more broadband pressures starting at lower frequencies (as low as 
200 Hz) than did the forward chord positions. Because of their broadband content, pres- 
sures from the aft chord positions did not provide good rotational- noise prediction results. 

Langley Research Center, 

National Aeronautics and Space Administration, 

Hampton, Va., October 9, 1974. 
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APPENDIX A 


DERIVATION OF WRIGHT’S SOLUTION 

This appendix presents a simplified derivation of the rotational noise equations 
developed in references 7 and 8 and outlined in the section "Rotor Rotational Noise 
Theory" of this, report. 


Generalized Solution of Wave Equation 

The rotor and force coordinate systems are shown in figure 20. The wave equation 
to be solved is 


1 9^p 9^p 

a^ 9t 2 9x 2 


9F i 

9Xj 


(x,t) 


(Al) 


where Fj[(x,t) = f j(x,t)6(x-y) is a general point force, y = y(t), the bar denotes a vector, 
and 5 ( ) is a Dirac delta function. 

The solution to equation (Al) then is (from Lighthill, ref. 30) 





9 

fi(y(t),t) 

dxi 

D 


where D = |x - y|, and £ J implies retarded time t - D/aQ. Therefore, 


p(x,t) 


- ^(x - y) 
47rD^ 


_ a 0 at D _ 


(A2) 


(A3) 


Solution of Wave Equation for a Rotating Point Force 
A rotating point force distribution can be modeled as follows: 

f(y(t),t) = f 07? (r)exp[j(7]f2t - W- 0)j (A4) 

where rj and £ are arbitrary phase constants. 
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APPENDIX A - Continued 


Hence, the radiated sound pressure for each rj is: 


S p - "ft" ^ + i'iexDM 
71 ' 4* D 2 M a 0 q ! ex P 3 


V®[t - - W - <p 

a o 


(A5) 


From figure 20, 


D = [r 2 + r^ - 2Rr cos a ' cos{0 - i /^| 2 


1 0rj 


Or] 




(-sin j3 sin ^), Xg ^ 11 /3 cos i^), Xg(-cos 0) 


where Xj, Xg* *3 are unit vectors in the x-p Xg, Xg directions, respectively. Now 
y • = 0 because fp acts in the plane normal to y. Therefore, 


(x - y) ■ f 0r] = x - = -R|f 0r? | [sin /3 cos o' sin(i p - 0) + sin a' cos /sj 


The sound radiation from a force system distributed over the rotor disk is deter- 
mined by letting j ^ 077 ( r )| = = P?]( r ) r dr d v^, which yields, for the radiation from 

the entire disk, 

SP^ = C f* P T7 (r) [sin /3 cos a' sin (\f/ - 9) + cos /3 sin o'] 

V J 0 J 0 47rD 2 11 


% + D) eXP < j 


0 


r dr di// 


(A6) 


The following simplifications are made in the analysis at this stage. In order to 
simplify the integration procedures, the observer azimuthal 6 is assumed to be zero 
since it will amount to only a phase shift in the final solution. It is assumed that the 
observer is in the far field. That is, D <-► R everywhere except in the phase term where 
D = R - r cos o' cos ij/. 

In addition, the following results are used: 
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z cos 


T)l//jjdl/' = 27Tj 77 J ?? (z) 



\jy- 7)tf/j]dJ// = - 2 277(-j) r? J T7 (z) 


Hence the sound pressure reduces to 


SP^ = A 


rP T j(r)^cos j3 sin o' sin /3^ 27r(-j)^J^ r cos cr'jdr (A7) 


?]S2r 


where |A| = - ~r~ exp 


4?rRa, 


0 


R 


* \ a 0 


jT70 t - 


exp(-j 0). 


If it is assumed that the loading is concentrated over Ar at an effective radius r e 
and Ar is small enough for to be constant over Ar, then 


SP. 


_ 7)12 
V 2Ra, 


0 


exp <] 


T)0[t - - (p j r e Pjj(re)Ar (cos /3 sin cr T 


sin J ^ Me cos 




(A8) 


The solution (A8) is for the sound radiation from a general pressure pattern, 
P 7? (r e )exp[j (r)fit - £\p - 0)1. The helicopter rotational-noise equations are developed along 
the same lines. All that remains to be done is to model the forces on the rotor disk. 
These forces are composed of uniform and fluctuating components. 

Uniform rotating forces .- For a helicopter the time history of the uniform rotating 
pressure pattern is as shown in figure 30. For ease of understanding, a rectangular dis- 
tribution over the chord at r = r e is shown, although any distribution can be used. This 
distribution can be represented as 

oo 

C^. exp(J(kBf2t - kBipjj (A9) 

k=-°° 
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APPENDIX A - Continued 


where 



¥§*k 


and 


X k = Chord spectrum function = 


sin kBc/2r e 
kBc/2r e 


for a rectangular distribution. 

Wright uses simple chord distribution functions such as rectangular, triangular, and 
half cosine in his theoretical development. For a particular sound harmonic m, the blade 
pressure mode of interest is 



Figure 30.- Uniform rectangular force f(t) 

in the rotor plane at position y; f(t) = C 

Cm = y • T = t' = — • t - 

. nr e x mB’ ft’ 0 ~ 


distribution seen by an observer on the fluid 

cos(mBflt - mBflt'); 
c 

nr e ’ 
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APPENDIX A - Continued 


P mB = -mB ex p[j( mBfit ' 


where P D = 

— mB 7rr e 


^mB anc * — mB * s ^dependent 0 f ^ f° r uniform forces. 
Fluctuating forces .- For fluctuating pressures, P g is dependent on \p, 


P mB = — mB^^ exp[j(mBfit - mB^/j] 


If the fluctuating pressures are constrained to have only periodic components, 


Z -m8,s exp(J(s^- 


s=-°° 


-mB + 


00 — 

1 -mB,s exp[ 5 (s^- + £ P mIM exp[j(s^- ^ 


s= 1 


s=-l 


where 


= — mB \ + ^ T ( ex P^ s ^“ 0s)| + expj^j (sip - 0 s j|J 


s=l 


p = BB - V 

— mB 7rr e *mB 


P mB,s —mB^s 


a s = 


= - ( P-< 


P = P 

- mB,-s _mB,s 


The pressure pattern for a particular sound harmonic mB and loading harmonic 


(A 10) 


that is, 


(All) 


s, is 
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represented as 


P mB = — mB ~2 




A 




where 


q ± = mB ± s (A13) 

Comparing equation (A12) with the general pressure pattern P^ gives 


r] = mB 

« = 


J_a JDCzc 

V r «> - 2^7 *«B 
0 = 0 S 


Hence, neglecting all phase terms, the sound radiation from a particular blade loading 
harmonic s is 


“s^mB 

SP„ n „ = -4 — — LB Ar 


— mB,s 2Ra, 


0 


q_ 


cos (3 sin o' — - ■ ■ sin j3 mBJ n (mBM e cos cr') 

mBM e q - 


“J. I 

+ cos /3 sin o' - „ sin (3 mBJ n , (mBM e cos o') 

mBMo l l+ e 


(A 14) 


Now, LB Ar = Lrp, Lp = L-p cos /3 is the total rotor lift force, and Dp = Lp sin /3 is 
the total rotor drag force. Hence, 

0!a 

SP = — v 

— mB,s 2 x rnB 


K L -K n ^V, + 


D mB/ r q. 


K L ' K D mBj y q 


(A15) 


N 


where K T == 

L Ra 


N 


Dt 


tt~> N = y a = mBJ a (mBM e cos o’), and 


J F Sin CT ’ D - Me’ " ~ 2i’ "q ■ 


JT g 

oig = — . Equation (A15) is the sound radiation from a particular loading harmonic s; it 
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does not include any phase terms. If all the phase terms are included, equation (A15) then 
becomes 



= SP 


mB,s 


exp 


jmBfi 



(A16) 


For convenience, 6 was equated to zero earlier in the analysis. As mentioned earlier, 

8 provides a phase shift of q0. For multiple blade harmonic radiation, that is, radiation 
due to all blade loading harmonics, the sound pressure at a particular sound harmonic mB 
is then 


SP mB = X SP mB,s 
s=0 


exp 




00 

1 


— mB,s 


(A17) 


J s=0 


Now, 



+ 



“s2±) y 

mB r q+ 


exp 



+ 



(A18) 


Substituting q for s as given by equation (A13) yields 
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zL — mB,s x mB X 2 1 Kl " mBj y q 6XP 


s=0 


q=mB 


K D q' 


( q0 r ^mB-qj ( -: > q+1 ) 


V VmB L K D q 

+ x mB Z “ T" 
q=mB 


K L - -SF/lq exp ' 


(« e + VmB)l (-i q + 1 ) <A19) 


Rearranging the terms in equation (A19) gives 


00 

1 — mB,s = T x mB( K L ' k d)%b exp^mBD- *„)_ 


• mB+l\ 

■J / 


s=0 


a. 


mB 


XmB( K L’'o) exp (' J VB) ( 'i ) 


I Ja T lSl *m B ( K L - + ( 2S q - ‘) «|q-n.Bj]J (-> q + 1 ) 


^ u |mB-q| 
q=-l 


K D q' 


2 x mB\ K L " m B/ y q 6Xp 


)( qe + ^|mB-q| 


-j q+1 


(A 20) 


where 


Sq = 0 for q < mB 


S a = 1 for q ^ mB 


The final result is obtained by substituting equation (A20) into equation (A17) after reversing 
the order of summation in the last term. 
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SP 


mB 


exp 


jmBfi 






mB0 - 



mB + 1 


+ "IT K L y 0 exp ('^raB) 


uu 

I 

q= l 


“|q-mB| 


K L- K D^)>- q ^P(i 


\i + ( 2S q ■ ') '’lq - niBl j 


* -^(Kl + K D ^)y q exp^qo- 4> |mB+q |)JJ-i q+1 ) 


(A21) 
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DETERMINATION OF OBSERVER AZIMUTH ANGLE 

The following experimental method was used to determine the blade- observer angle 
6. During data acquisition, a blade- position indicator was recorded on tape. The indica- 
tor consisted of a large pulse followed by three smaller pulses (fig. 31(a)). The large 
pulse indicated when the instrumented blade passed a fixed point on the rotor disk plane. 
The smaller pulses divided the blade path into quadrants. 

Let point A (fig. 31(b)) be a fixed point on the disk plane and let M be the micro- 
phone (observer) position. Then, 9' is the angle between lines OA and OB corresponding 
to t' in the rotor plane, and <p is the angle between lines OA and OM. Point B repre- 
sents the time at which data sampling for Fourier analysis was stated. Then, the value of 
8, the blade- observer angle at the onset of data analysis, is 

9 = 4> - e' 


Data channel 


Position indicator 
► t 

(a) Sample oscillograph record, (t’ is defined in eq. 31(b).) 
Figure 31.- Determination of 9, the blade- observer azimuthal angle. 
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M Microphone position 

A Blade position marker 

B Starting point for data 

t' Time taken for the blade 
to travel an angle 6 ' 
from A to B 



M 


(b) Observer orientation. 
Figure 31.- Concluded. 
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